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graph: an explicitly solvable model 

Stanislav Molchanov and Boris Vainberg 



Abstract. We consider an explicitly solvable model (formulated in the Rie- 
mannian geometry terms) for a stationary wave process in a specific thin do- 
main Q E with the Dirichlct boundary conditions on dQ E . The transition from 
the solutions of the scattering problem on Q E to the solutions of a problem on 
the limiting quantum graph T is studied. We calculate the Lagrangian gluing 
conditions at vertices v £ V for the problem on the limiting graph. If the 
frequency of the incident wave is above the bottom of the absolutely continu- 
ous spectrum, the gluing conditions are formulated in terms of the scattering 
data of a problem in a neighborhood of each vertex v 6 T. Near the bottom of 
the absolutely continuous spectrum the wave propagation is generically sup- 
pressed, and the gluing condition is degenerate (any solution of the limiting 
problem is zero at each vertex). 
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1. Introduction 

The paper concerns the asymptotic analysis of the wave propagation through a 
system of wave guides (fibers) when the thickness e of the wave guides is very small 
and the wave length is comparable to e. The simplest model one can consider is 
the stationary wave (Helmholtz) equation 

(1.1) -Au = -4u, x G Q £ , 

in the domain Sl e C R v 1 v > 2, consisting of finitely many cylinders (tubes) of 
lengths li, £2, • • - In with the diameters of the cross-sections of order O (e). Some of 
the lengths can be infinite. The junctions Ji, • • • , Jjv' connecting the cylinders into 
networks are compact domains in R" with the diameters of the same order 0{e). 
The axes of the cylinders and the centers of the junctions form edges and vertices, 
respectively, of the limiting (e — > 0) metric graph T. 

The Helmholtz equation in Vt e must be complemented by the boundary condi- 
tions (BC) on dVl e . In some cases (for instance, in a study of heat transport in VL e ) 
the Neumann BC is natural. In fact, the Neumann BC presents the simplest case 
due to the existence of a simple ground state (a constant) of the problem in f2 e . 
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However, in many applications, the Dirichlet, Robin or impedance BC are more 
important. We will consider (apart from a general discussion) only the Dirichlet 
BC, but all the arguments and results can be modified to be applied to the problem 
with arbitrary BC. 

We are going to study an explicitly solvable model. We assume that all the 
tubes have the same cross sections ui e . We also assume that oj e is e— homothety of 
a domain w S R v ~ l . Let Ao be the principal eigenvalue of the Laplacian Hq = 
— A„_i in lu. Thus, e~ 2 Ao is the principal eigenvalue of — A„_i in uj e . In the 
presence of infinite fibers, the spectrum of the Dirichlet Laplacian on Q £ has an 
absolutely continuous component which coincides with the semi-bounded interval 
[e~ 2 Ao,oo). The equation (| 1 - 1|> is considered under the assumption that A > Ao, 
when propagation of waves is possible. There are two very different cases: A — * Ao as 
e — > 0, i.e. the frequency is at the edge (or bottom) of the absolutely continuous 
spectrum, or A — > A > Ao, i.e. the frequency is above the bottom of the absolutely 
continuous spectrum. 

If e — ► 0, one can expect that the solution u e of (jl.l|l on fl e is close to the 
solution u = u e of a much simpler problem on the graph T. The function u satisfies 
the following equation on each edge of the graph 



where s is the length parameter on the edges. One has to add appropriate gluing 
conditions (GC) on the vertices v of T. These gluing conditions give basic informa- 
tion on the propagation of waves through the junctions. They define the solution 
u of the problem (|1.2f> on the limiting graph. The ordinary differential equation 
(11.21) . the GC, and the solution u depend on e. However, we shall often call the 
corresponding problem on the graph the limiting problem, since it enables one to 
find the main term of small e asymptotics for the solution u = u £ of the problem 

oil in n e . 

The convergence of the spectrum of the problem in fl e to the spectrum of 
a problem on the limiting graph has been extensively discussed in physical and 
mathematical literature (e.g., pQ-jl], Q5J 1101 1121 II 3j and references therein). 
This list, containing important contributions to the topic and some review papers, 
is far from complete. What makes our paper different is the following: all the 
publications that we are aware of, are devoted to the convergence of the spectra 
(or resolvents) only in a small (in fact, shrinking with s — ► 0) neighborhood of 
Ao (bottom of the absolutely continuous spectrum), or below Ao- Usually, the 
Neumann BC is assumed. We deal with asymptotic behavior of solutions of the 
scattering problem in Sl e when A is close to A > Ao, and the BC on <9f2 e can be 
arbitrary. It turns out that the GC on the limiting graph in the case A — ► A > Ao 
is different from the case when A — ► Ao. 

Papers [2], |1U| . |13) contain the gluing conditions and the justification of 
the limiting procedure e — ► in the case when the Neumann BC is imposed at the 
boundary of fl e , and A — > (Ao = in the case of the Neumann BC). Typically, the 
GC in this case are: the continuity of u (s) at each vertex v and Y]j—i u'j (v) — 0, i.e. 
the continuity of both the field and the flow. These GC are called Kirchhoff 's GC. 
In the case when the shrinkage rate of the volume of the junction neighborhoods 



(1.2) 
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is lower than the one of the area of the cross-sections of the guides, more complex 
energy dependent or decoupling condition can arise (see [10] . [4] for details). 

Let us stress again that this is the situation near the bottom Ao = of the 
absolutely continuous spectrum. A recent paper by O. Post |12| contains analysis 
of the Dirichlet Laplacian near the bottom of the absolutely continuous spectrum 
Ao > under the condition that the junction is more narrow than the tubes. It was 
proved in |12j that in this case the GC for the problem on the limiting graph are 
the Dirichlet conditions, i.e. waves do not propagate through the narrow junction 
when A is close to the bottom of the absolutely continuous spectrum. 

Our paper concerns the asymptotic analysis of the Dirichlet problem when 
A is close to A > Ao. As we shall see, the limiting GC conditions in this case 
differ from the Dirichlet or Kirchhoff's conditions. They are formulated in terms 
of the scattering coefficients of the original problem. At the first glance, it looks 
like a vicious circle: one needs to solve the problem in il e and find the scattering 
coefficients in order to state the limiting problem. However, this is not the case. 
The scattering coefficients arc needed only for local problems in domains which 
consist of one junction and adjoint tubes, but not for the problem in f2 e . 

We emphasize again that we do not consider equation in a branching do- 
main in Euclidean space, but rather a similar equation on a product of a compact 
manifold and a quantum graph. The latter equation admits a separation of vari- 
ables. The proposed model does not have direct physical significance. However, the 
changes in the metric (or potentials) in our model play the role of junctions, and 
the problem under consideration is a simplified analog of a more realistic situation. 
The present model is chosen only for the sake of simplicity. We proved that the 
obtained results are valid for the problem in more general domains fl e (where the 
separation of variables is impossible). This general case will be discussed elsewhere. 

We also studied the Dirichlet problem for general domains fi e without spe- 
cial assumptions on the geometry of the junctions when, simultaneously, e — ► 0, 
A — > Ao, and the diameters of the guides and junctions have the same order 0(e). 
Our conclusion is that, generically, the limiting GC in this case is the Dirichlet 
condition. Thus (generically) waves will not propagate through the junctions when 
the frequency is close to the bottom of the absolutely continuous spectrum. Let 
us stress that this is true both in the case when the diameters of the junctions are 
smaller than the diameters of the guides, and in the case when they are larger. 
Some special conditions (they will be described in another paper) must be satisfied 
for waves to propagate if A — > Ao. An infinite cylinder, which can be considered as 
two half-infinite tubes with the junction of the same shape, can be considered as an 
example of a domain where the propagation of the wave at A = Ao is not suppressed. 
Practically, we do not deal with the problem near the bottom of the absolutely con- 
tinuous spectrum in this publication. A detailed analysis of this problem will be 
published elsewhere. However, we show here that the GC on the limiting graph 
for our simplified model with A = A, generically, have a limit as A — > Ao and the 
limiting conditions are the Dirichlet conditions. 

The authors are grateful to P. Exner, P. Kuchment, and S. Novikov for very 
useful discussions. 
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2. Setup of the model. 

To understand better the general model under consideration we start with a 
simple example. Let us consider the Helmholtz equation 

(2.1) -Au= 4", u\ dn = 

inside of a thin 2-dimcnsional strip f2 e which has a constant width e outside of the 
e-neighborhood of the origin. This strip is a particular case of domains f2 e described 
in the introduction. In our case Q £ consists of two half strips with one junction. 
Consider the case e = 1. We roll up the strip £7 = fl\ and transform it into a 
cylindrical surface C with circle cross-sections. The radius of the cross sections is 
constant outside of the neighborhood of the origin. Let (ip, s), < ip < 2tt, s G i?, 
be cylindrical coordinates on C, and let A(s) be the radius of the cross-section of 
the cylinder C. The boundary d^l corresponds to a cut of the cylinder along the 
meridian ip = 0. 

We replace the Euclidean metric on C by the metric of the surface of revolution: 
dl 2 = ds 2 +A 2 (s) dip 2 , and we replace the Laplacian A on by the Laplace-Beltrami 
operator on C which is defined by the metric dl 2 . Then we arrive at 

d 2 u A' (s) du 1 d 2 u , . . . . 

W ~W - aU^ - A^W = XU > U(M = U(2?r ' S) = °' 
For arbitrary e > 0, we assume that il e is e-homothety of fi, i.e. f2 e is the image of 
Q e under the action of the map (s,ip) — > e~ 1 (s,ip). Then the problem (|2.2|l takes 
the form 

d 2 u A' (s/e) du 1 d 2 u X . . . . . 

(2 ' 3) -^-7Xm^s-^mW = ^ ^s) = u^/e,s)=0. 

One can expect that the propagation of waves governed by this modified equation 
is very similar to the one governed by the original model (|2.1|1 . These two models 
also have similar levels of difficulty. 

Our general model is a direct generalization of (|2.2() . This model is still very 
special and does not include most of branching domains f2 e in Euclidean spaces. 
However, we believe that the model captures some qualitative characteristics of 
wave propagation in thin structures. The main advantage of the model is that it 
allows the separation of variables. 

Let F be a connected quantum graph (see with a finite number of 

vertices {u,} = V and edges {ej} = E. It is assumed that each edge ej has a 
parameterization s € (0,lj), 1 < j < N, and some edges can have infinite length. 

d 2 

One can introduce the basic Hamiltonian Hq = — — on the space C^°(r\y) 

ds 2 

of infinitely smooth functions, supported outside the set V. The space L 2 (T) of 
square integrable functions on T and the Sobolev space H 1 (r) can be introduced 
in a natural way. Later we will introduce self adjoint extensions of the operator Hq 
which are defined by GC at vertices v e V. 

Let M be a Riemannian (y — l)-dimensional manifold with the Riemannian 
metric dip 2 and nontrivial boundary dM . Let M e be the image of M under the 
action of the map <p —> <p/e. The equation will be given on the manifold il e = T x M £ 
which is a Cartesian product of a quantum graph T and M £ . We change the metric 
on Vt s in order to take into account junctions between tubes ej x M e . It is more 
convenient to introduce first the metric on the rescaled manifold £1 = T £ x M where 



NETWORK OF THIN FIBERS 



5 



edges e £ of V s have lengths lj/e. The Riemannian metric dl 2 on f2 has the same 
form as in the example above, i.e. for each edge £ P, 

(2.4) dl 2 = ds 2 + Aj (s) 2 dip 2 , s G (0, lj/e) 

Here Aj are smooth positive functions on e £ such that Ay(£) = 1 for £ G (1, — 
1). The Laplace-Beltrami operator on associated with the metric dl 2 has the 
following form on each part x M of f2 : 

. a 2 ^ ( s ) a i . 

A = I J I A 

ds 2 Aj (s) ds A? (s) Vl 

where A v is the Laplace-Beltrami operator on M. Thus, the equation —Am = Xu 
on Q, after the inverse rescaling, corresponds to the following equation on £! e : 

d 2 u Aj (s/e) du 1 A 

(2 ' 5j eAjia/e)^ A^Jsje) U ~ s 2 ™' &j E 

One could arrive to the equation (|2.5|l by introducing the appropriate Riemannian 
metric directly on f2 e without rescalings made above. 

Equation i|2.5f) is symmetric with respect to the measure da = Aj(sj 'e)dsdip. 
The equation can be simplified (reduced to an equation, which is symmetric with 
respect to the measure dsdip) using the substitution 

u(s, if) — ► Aj (s/e)u(s,if), e 3 eE. 

The new equation (we preserve the same notation u for the solutions after the 
substitution) has the form 
(2.6) 



d 2 u 1 



8s 2 e 2 A 2 (s/e) 



A'!(s/e) lA'j(s/e)^ 



E 



\u, ej e E. 



2A 3 {s/e) ^Aj{s/e)' 

The equations l|2. 6f) have to be complemented by the appropriate BC, GC and 
conditions at infinity. We impose the Dirichlct BC on T £ x dM : 

(2.7) u = on T e x dM. 

Let us recall that the degree d = d(v) of a vertex v is the number of edges with an 
end point at v. We split the set V of vertices in two subsets V = V\ U V2, where 
the vertices from the set V\ have degree 1 and the vertices from the set V2 have 
degree at least two. We impose an arbitrary symmetric homogeneous BC on v x M, 
v G Vi, for example the Dirichlet BC: 

(2.8) u = on v x M, v G V\. 

One has to specify the meaning of the equation l|2.6|) on v x M . We assume 
that the solution u of (|2.5|) satisfies Kirchhoff's GC at vertices v G V2 , i-e. u is 
continuous at each vertex v G V2 and the flow is preserved: 

d(v) „ 

(2.9) «6C(fl E ), E^( u ) =0 > weF2 - 

Here is the differentiation with respect to the parameter s on the edge ej in the 
direction out of the vertex v. The conditions l|2. 9f) arise often in applications. The 
conditions l|2.9|) are also in agreement with the GC in our first example, where fl e 
is a strip and the functions on the strip are smooth. Let us stress that Kirchhoff's 
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GC at vertices v € Vi hold for the solution of the problem in f2 e . As we will see, the 
main term of asymptotics of that solution, which is determined by the equations 
on r, satisfies different GC. 

In order to complete the statement of the problem (|2.6fl - l|2.9ll . one needs to 
specify conditions at infinity. We are going to consider eigenfunctions and solutions 
of the scattering problems. To describe the latter solutions, we note that the 
problem l|2.5|l - l|2.9[l admits the separation of variables: 

(2.10) u( 7 , ip) = 1>(>y)a(<p/e), 

where 7 is a point on the graph T and a(ip) is an eigenfunction of the operator A v : 

-A v a((p) = X'a(tp), <p e M; a{ip) = on dM. 

Obviously, a(tp/s) is an eigenfunction of the Dirichlct Laplacian on M E with the 
eigenvalue A'/e 2 . Then tp satisfies the equation 

(2.11) + s- 2 Q(s/£)ip= (^) 2 V>, fc= v / A^A 7 >0, 
where 

A" 1 A' 1 — A 2 

Q = Qi = 2A 3 - i ( i; )2 + A ' (_ 4T ) on 6j e E ' 

and the boundary conditions 

(2.12) V = at «6Vi, 

(2.13) V G C(T), ^^( u)=0 , veV 2 . 

j = l 3 

Let E' C E be the set of semi-bounded (lj — 00) edges {ej}, 1 < j < r. We 
have a natural parameterization s £ (0, 00) on ej € i?' where s — corresponds to 
the end point of ej. Let us recall that the coefficients Aj differ from 1 only in 
neighborhoods of the vertices, i.e. Aj(s) = 1, s > 1 if j < r. Then the equations 
(|2.11|l imply 

-ip" = -^ip on e, e E', s > e; fc = VA - A' > 0, 

i. e. V = ae jfcs / £ + f3e~ lks / £ on ej S when s > e. 

Definition 2.1. The function V^ m \ 1 < m < r, is called a solution of the 
scattering problem (on the graph) if it satisfies l|2.11|l - H2.13ll and 
(2-14) 

M = ( e-^+J^^ ° n e ™ EE '> .* >£ , k = VA^ > 0. 
t Tj.mt - 1 on ej € E, s > e, j ^ m 

The function u = given by p.lOjl with tp — ijj^ is called the solution of the 
scattering problem in Q £ . These functions describe the propagation of a wave of 
unit amplitude incoming through the wave guide e m x M, r m m is the reflection 
coefficient, Tj jm , j ^ m, are the transmission coefficients. 

Due to the separation of the variables, the scattering problem in Q E is reduced 
to the scattering problem on the graph. However, this problem on the graph is 
governed by the equation (|2.11|l with variable coefficients. Obviously, the potential 
e~ 2 Q(s/e) is supported in e-neighborhoods of the vertices. 
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Finally, wc come to the main object of the investigation: asymptotic analysis 
as e — > of the scattering solutions u — w- m \ In fact, the next section contains 
asymptotic analysis of arbitrary solutions of the problem H2.6|H2.8fl f° r which sep- 
aration of variables i|2.10|l holds. For example, u can be an eigenfunction of the 
problem. Obviously, formula (|2.10fl reduces the analysis of the asymptotic behavior 
of the function u to the study of the asymptotic behaviour of the solutions ip of the 
problem 

3. Asymptotic behavior of solutions. 

Let ip be an arbitrary solution of the problem (|2.11() - (|2.13() . The main feature 
of the Schrodinger equation Q2.11|l is that the potential is very singular (of order 
e~ 2 ) near the vertices and it vanishes outside of a very narrow (of the size e) 
neighborhood of the vertices. Then ip has the form ip = otje lks / £ + (3je~ lks l £ on each 
edge ej € E outside of e-neighborhoods of the vertices. We will call the function 
ip on r, which is equal to aje lks t £ + /3je~ lks ' e on each edge ej € E up to the 
end points, the "limiting" function. Let us stress that the "limiting" function still 
depends on e, and this is the reason to use the word "limiting" in quotation marks. 
The "limiting" function provides a good simple approximation for the solution ip 
since it differs from ip only in e-neighborhood of the vertices. The goal of this section 
is to find a way to determine the "limiting" functions for scattering solutions and 
for eigenfunctions of the problem on the graph directly (without solving the original 
problem). 

The "limiting" function ip for any solution ip of l|2.11|l ) satisfies the equation 

_V/' = (-) 2 V; onT. 

e 

However, if ip satisfies the GC (|2.12|) . the GC for ip are different from l|2.12|l . 
In order to find the GC for ip one needs to study supplementary problems for 
simple subgraphs T v of T, consisting of one vertex v of degree d = d{v) and the 
edges ej 1 ,...,ej i with an end point at this vertex. Let us consider one of such 
supplementary problems, when d(v) > 1 (i.e. when v G V%). 

We fix the parameterization on the edges e^, e Jd in such a way that s > 
and s — corresponds to the point v. Recall that the potential e~ 2 Q(s/e), defined 
on L, vanishes outside e-neighborhoods of the vertices. Let the potential e~ 2 B(s/e) 
on coincide with e~ 2 Q(s/e) in the e-neighborhood of the vertex v and B = 
when s > e. Let (p( m > be solutions of the scattering problems for the Schrodinger 
operator on the subgraph T v with the potential B : 

(3.1) -ip" + e- 2 B{s/e)ip = (-)V on r„, 

£ 

(3.2) <p e c(r„), E^W = o, 

(3.3) ip = p^ = 8 n , m e~ lks / £ + t n , m e lks ' £ on e jn , s>£. 

Here tp n is the restriction of p on ej n , S niTn is Kronecker's symbol. Of course, 
coefficients t n>m depend on v. Let 

Tv — [^n,m] 
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be the matrix of scattering coefficients. The diagonal elements t m ,m of the matrix 
T v are reflection coefficients of the corresponding scattering solutions tp 1 --" 1 ^ , and the 
other elements are transmission coefficients for <p( m ) . 

For any solution (p of l|3.1l) and its "limiting" solution 92, we denote by <p the 
column vector whose coordinates are restrictions tp n of tp on ej n . 

Theorem 3.1. Let d(v) > 1 (i.e. v S V 2 ). Then 

1) The matrix T v — [t n _ rn ] does not depend on e. 

2) For any function ip which satisfies the equation 13. 1\) and Kirchhoff's' GC 
\3.ty . the "limiting" function tp satisfies the following GC at the vertex v: 

(3.4) l J. {I + Tv )^-(I-T v )p = 0, 

tz as 

where I is d x d unit matrix. 

3) The matrix T v is unitary and symmetric (t n ^ rn = t m . n ); the d x 2d matrix 
+ T v ), (I — T v )] has rank d. 

Remarks. 1) If the matrix I + T V is not degenerate, then i|3.4|) can be written 
in the form 

f = Cp, C=—(I + T v y\l-T v ), 
ds e 

Here the matrix C is Hermitian (C = C*) due to the unitarity of T v . If / — T v is 
not degenerate, then tp = C -1 ^ where C" 1 is Hermitian. 

2) The scattering matrix T v can be expressed through the scattering solutions 
of simpler scattering problems on individual edges, see below. 

3) The unitarity of the scattering matrix is a standard fact in the scattering 
theory. Its symmetry is also well known for 1-D Schrodinger equation. In 1-D case 
it means that the transmission coefficient does not depend on the direction of the 
incident wave. The authors learned about the symmetry of the scattering matrix 
in a more general situation from S. Novikov (see 

Proof. In contrast to the case of the whole graph T, e-independence of the 
scattering matrix T v for a star shaped graph r„ becomes obvious after the rescaling 
s — * se. 

Let t^™ 1 ' be the " limiting" function which corresponds to the scattering solution 
ip( m \ i.e. 

<? (m) = S n>m e- iks/e + t n:m e lks/e on e jn for all s > 0. 
Then the Cauchy data for the corresponding vector </}( m ) at the vertex v are 

fi (m) (v) = S m + t m , ^j—{v) = —(-6 m + tm), 

ds e 

where 6 m is th.G column vGctor with coordincites 5 n rni 1 ^ ft <C and t m is the 
m-th column of the matrix T. Hence, for the matrix $ with the columns ^ m ^(s), 
1 < m < d, we have 

$(v) = I + T v , ^( v ) = -(-I + T v ). 
ds e 

This immediately implies that the "limiting" functions of the scattering solutions 
<^ m ) satisfy the GC 

It is obvious that the solution space for the problem l|3.1Jl , (|3.2|l is d-dimensional 
and that the scattering solutions are linearly independent. Thus any solution ip of 
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(|3.1(1 . 1(3. 2|l is a linear combination of the scattering solutions, and therefore its 
"limiting" function <p satisfies the GC ((3.4(1 . 

Let us prove the last statement of the theorem. If some of the edges of the 
graph T v are finite, one can extend them to infinity. Thus, without loss of the 
generality we can assume that all these edges are infinite. Let Ty be the part of 
r„ on which s < a. Consider two scattering solutions (p( mi >, (p( m2 > on T v defined 
by 1(3. 1(1 - 13. 3|) . Since the problem ((3.1(1 . 1(3.2(1 is symmetric, from Green's formula 
for Lp {m ^ and Tp^ on T ( v a) it follows that 

d 

E 

n=l 

If we substitute 1(3. 3|) in the last formula, we arrive at 

d 

Ef 7 _ 7 -lika/e , . 2ika/e _ r _ r> 

ri=l 

We take the average with respect to a £ (A, 2 A) and pass to the limit as A — ■> oo. 
This leads to the orthogonality of the columns t mi , t„ l2 of the matrix T v when 
mi 7^ TO2 and to the condition \t mi \ = 1 when mi = rri2- Thus, the matrix T v is 
unitary. The symmetry of T v can be proved similarly using Green formula for ip( mi ^ 
and ip( m2 \ Since the second part of the third statement of the theorem is obvious, 
the proof of the theorem is complete. □ 

Let now v G V% and let e v be the edge of T with an end point at v. One 
can assume that the parameterization on e v is chosen in such a way that s = 
corresponds to the vertex v. 

We make the rescaling s — > se on e v and consider the following supplementary 
problem on the rescaled edge: 

-(p" + Q(s)(p = k?tp, s<--l; ¥>(0)=0, ^(0) = 1, 

£ 

where = \e v \ is the length of e v . Then 

<p = a v e~ lks + p v e lks , s £ (1, — - 1), 

e 

where a v , (3 V do not depend on e. 

Let ifj satisfy H2.11|l on e v and tp(v) = 0. Then ip(s) on e«, s < l v — e, is 
proportional to ip(s/e). Hence, the "limiting" function -0 has the following form on 

e v , v G Vl : 

i> = C{a v e~ lks/£ + P v e lks/e ) on e v , s < l v - e, 

and therefore 

%£ — — 

(3.5) -(a v +f3 v )ij'(v)-(a v -(3 v )ip(v)=Q ) v G V x . 

We note that Theorem 13.11 and the formula l|3.5|l provide a local description 
of the GC for "limiting" functions, where the solution can be defined only in a 
neighborhood of individual vertices. Thus, the following theorem is an immediate 
consequence of the Theorem 13.11 and the formula 1)3.5(1 . 



7 (mi) 

as 



(m 2 ) 



(a) = 0. 
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Theorem 3.2. Let the function ip satisfy H2.ll)) and the Kirchhoff's' GC 12.1 
Then the corresponding "limiting" function ip is a solution of the equation 

(3.6) = (-) 2 ^ on r, 

e 

which satisfies $3. 5)) and the GC 

(3.7) l JL( I + Tv )^-(I-T v )$=0, veV 2 . 

In particular, if ip is a solution of the scattering problem 12.11)) - 12.13)) . then ip 
is the solution of the scattering problem for the equation $3. 6)) with the conditions 
$3.5)) and |ff.7| ). If ip * s an eigenfunction of the problem $2.11)) . $2.12)) . then ip is 
an eigenfunction of the problem \S.f>)) . \3.5)) . \3.T[ . 

4. Evaluation of the scattering matrix T v . 

In the last part of the paper we will reduce the scattering problem (|3.1|l - (|3.3|) 
on the star shaped graph T v to d simpler scattering problems on individual edges 
ej n extended to infinity in both directions. Let the straight line R with coordinates 
s G (—00,00) represent the extended edges ej n . Let the potential B n (s) on R 
coincide with the restriction of B(s) on ej n when s > 0, and B n (s) = for s < 0. 
Let ip = %p( n ) , 1 < n < d, be the solution of the following scattering problem on the 
line R (the extension of the edge ej n ) : 

-ip" + B n (s)ip = k 2 ip, - 00 < s < 00; 
(4.1) V = e iks + r n e~ lks , s < 0; ip = t n e iks , s > 1. 

This problem describes the propagation of the wave incoming from s = — 00. 

The following theorem allows one to express the scattering data t n _ rn of the 
problem (|3.1I) - (|3.3|) on the graph T v through the reflection r n and transmission t n 
coefficients of the problem (14. IB . Let 

d 

1 - r-j 



f-f 1 + r-j 

J= i 

Theorem 4.1. The following formulas hold: 

(AO) t - 2tmtn 777 ^77- t - ^ tjnQ + Tm) 

(l + r m )(l + r n )p (l + r m yp t m (l + r m ) 

Proof. Let us look for the solution <p( m ) I3.1|) - (|3.3|l on the graph in the form 

fA*\ »/ tl _f 7n^(s/e) one,„, n^m 

( dj * ~ I ^ m {s/e)+ lm iP m (s/e) one jm 

with constants 7„ which will be chosen below. In fact, the constants "f n = 7n,m 

de- 
pend also on to, but we will often omit the index to to simplify formulas. Obviously, 
<^ m ) satisfies (JXTJ and (ESJ with 

(4.4) t nm — ^in.mXn- 

It remains only to choose j n = 7n !m in such a way that i|3.2|) holds. Then (|4.3|l will 
be the scattering solution of the problem on the graph. 
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From 14. 1|) it follows that 

4> n (0) = l + r n , =*(-l + r n ). 

Thus, the substitution of 1)4.3)) into l|3.2() leads to the following system 
(4.5) 

1 + f 1 r d 

7n(l + r n ) = j ~ + 7m (1 + r m), n ^ 771 ; — + ^ 7n(~l + T*n) = 0- 

£m ■ 

Hence 

_s l+f m 7m0- + r m ) 

(4 ' 6) 7 " = 7 "^UlT^ + (i + rn ) ' n + m - 

We substitute this expression for j n into the last equation 1)4.5)1 and arrive at the 
equation for 7 m : 

l-f m ^ {l+f m ){-l + r n ) 



(l + r m )p 7m =0, 



which is equivalent to 

, . „v l-r m (l + r m )(l-r m ) 1 + F m , . 

(4-7) + = p - (1 + r m )p7 m = 0. 

The sum of the first two terms is equal to 
(l-r m )(l + r m ) + (l + r m )(l-r m ) =2 _1- |r m | 2 _ 2\t m \ 2 _ 2t r , 



t m {l + r m ) t m (l + r m ) t m (l + r m ) (l + r m ) 

From here and l|4.7|) it follows that. 

_ 2£m 1 -|- T m 

m ' m (l + r m ) 2 p t m (l + r m )' 

This together with 1)4.6)) and (|4.4|) implies 1)4.2)) . The proof is complete. □ 

5. Solutions near the bottom of the absolutely continuous spectrum. 

Recall that is the bounded part of the star shaped graph T v on which 
s < a. When k = \/X — A' > 0, the GC at vertices v £ T for the limiting problem 
on the graph are formulated (see Theorem l3.1|) in terms of solutions of the scattering 
problem (|3.1I) - (|3.3() . In order to describe the behavior of the solutions when k = 
VA — A' — * 0, we need to consider the Neumann problem on with the BC at 
points where s = a : 

(5.1) -ip" + e- 2 B(s/e) V = Q onrW, 



(5.2) peCCTM), E^H=0; ^=0, l<n<d. 

3 = 1 

The choice of a above is not important since B(s/e) = for s > e and solutions 
of 1)5. 1)5. 2)1 are constants when s > e. For example, one can choose a = 2/e. 
One can also simplify the problem 1)5.1(1 . 1)5.2)1 by changing the variable s — > se. 
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The function ip in the new variable is an eigenfunction with zero eigenvalue of the 
following Neumann problem: 

(5.3) -if" + B(s)<p = onlf), 



(5.4) ?€C(TW), J2^(v) = 0; = 0, l<n<d. 

Zero is not an eigenvalue of the problem (|5.3() . 1)5.4)) for a generic potential 
B(s). Thus, the conclusion of the next theorem holds generically. 

Theorem 5.1. Let zero be not an eigenvalue of the problem \5.$ - Then 

the GC \3-4\l respectively) has the limit as k — ► 0, and the limit condition is 

the Dirichlet condition: <p n (v) — 0, 1 < n < d. 

Proof. Consider the restriction of the scattering solution ip = p^ of the 
problem ()3. 1)1 - 1)3.3)1 to t{?\ a = 2/e. After the rescaling s — * se, the following 
relations hold 

(5.5) -ip" + B(s)cp = k 2 ip onr< 2) , 
(5.6) 

p€C(T< >), J2^( v )=0 ; ^^ = -i*(<y n , m -t nim ) ) l<n<d. 

3 = 1 

(2) 

Since |t n , m | < 1, from these relations it follows that ip — O(k) onT„ as k — > 0. In 
particular, </? n (2) = 0(k) as k — > 0. The latter together with 13 . 3|) implies that 



Differentiation of 1)5. 5[l . 1)5.6)1 with respect to k allows one to conclude that ^jff-(2) = 
0(k) as k — > 0. From here it follows that Juin.m = 0{k) as fc — * 0. Hence, I — T = 
0(k), JrT = O(fc), as jfe -> 0. This allows one to pass to the limit as k — > in 1)3.4)1 . 
()3.7I) and get the Dirichlet GC for the limiting function. The proof is complete. □ 
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